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ABSTRACT 

Damped Lya Systems (DLAs), with N(HI) > 2 x 10^" cm^^, observed in the 
spectra of quasars have allowed us to quantify the chemical content of the Universe 
over cosmological scales. Such studies can be extended to lower column densities, in 
the sub-DLA range (10"^^ < N(HI) < 2 xlO^*^ cm~^), which are systems believed to 
contain a large fraction of the neutral hydrogen at z > 3.5. In this paper, we use a 
homogeneous sample of sub-DLAs from the ESO UVES archives presented in Paper I 
(Dessauges-Zavadsky et al. 2003), to observationally determine for the first time the 
shape of the column density distribution, f{N), down to N(HI) = 10^^ cm~^. The 
results are in good agreement with the predictions from Peroux et al. (2003). We 
also present the kinematic and clustering properties of this survey of sub-DLAs which 
appear to be marginally different from the DLAs. We compare low- and high- ionization 
transition widths and find that the sub-DLAs properties roughly span the parameter 
space of DLAs. We also find hints of an increase of metallicity in systems with larger 
velocity widths in the metal lines, although the statistical significance of this result is 
low. 

We then analyse the chemical content of this sample in conjunction with a compi- 
lation of abundances from 72 DLAs taken from the literature. As previously reported, 
the individual metallicities traced by [Fe/H] of these systems evolve mildly with red- 
shift. Moreover, we analyse the HI column density-weighted mean abundance which is 
believed to be an indicator of the Universe's metallicity. Although the number statis- 
tics is limited in the current sample, the results suggest a slightly stronger evolution 
of this quantity in the sub-DLA range. The effect is predominant at z < 2 and most 
of the evolution observed lies in this redshift range. Observational arguments support 
the hypothesis that the evolution we probe in the sub-DLA range is not due to their 
lower dust content. Therefore these systems might be associated with a different class 
of objects which better trace the overall chemical evolution of the Universe. Finally, 
we present abundance ratios of [Si/Fe], [0/Fe], [C/Fe] and [Al/Fe] for sub-DLAs in 
conjunction with DLA measurements from the literature. The elemental ratios in sub- 
DLAs are comparable with those from DLAs. It is difficult to decipher whether the 
observed values are the effect of nucleosynthesis or are due to differential dust de- 
pletion. The metallicities are compared with two different sets of models of galaxy 
evolution in order to provide constraints on the morphology of quasar absorbers. 
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1 INTRODUCTION 

In addition to traditional emission studies, absorption sys- 
tems along the line-of-sight to distant quasars provide a com- 
pletely independent probe of galaxy evolution. The highest 
column density Damped Lya systems (hereafter DLAs) have 
N(HI) > 2 X 10^° atoms cm . The reason why damped sys- 
tems are a cosmologically important population is that they 
contain most of the neutral gas in the Universe at « > 1 
(Lanzctta et al. 1991; Wolfe ot al. 1995; Storric-Lombardi 
et al. 1996b). Since the metal content of these systems can 
be determined with rather high precision up to high red- 
shift, they are a powerful tool to study the chemical evo- 
lution of galaxies. In particular, a way to trace the metal- 
licity of the Universe is provided by estimating the ratio of 
the total metal content to the total gas content measured 
in these systems (Pcttiui ct al. 1997; Pel, Fall & Hauser 
1999). Using such techniques, Prochaska & Wolfe (2002), in 
agreement with previous work, find that the hydrogen col- 
umn density- weighted [Fc/H] mctallicitics have similar val- 
ues from 1.5 < z < 3.5. The results arc in disagreement with 
expectations since models of cosmic chemical evolution (e.g. 
Pei & Fall 1995) predict evolution of global metallicity with 
cosmic time. Savaglio (2000), however, has used 50 DLAs 
and lower column density systems and finds clear evidence 
for redshift evolution, although this sample suffers high in- 
homogeneity. 

In a recent study, Peroux ct al. (2001) have shown that 
sub-DLAs, defined as systems with 10^" < N(HI) < 2 x 10^° 
atoms cm~^, play a major role, especially at higher redshifts. 
The authors postulate that at z > 3.5, 45% of the neu- 
tral gas mass is in sub-DLAs. These predictions are based 
on a constrained extrapolation of the quasar absorber col- 
umn density distribution, f{N), to lower column densities 
assuming the distribution can be fitted by a gamma func- 
tion. This suggests that the metallicities of sub-DLAs re- 
quire detailed study in order to obtain a complete picture of 
the redshift evolution of the metallicity of the Universe. For 
this purpose, we have constructed a homogeneous sample of 
sub-DLAs based on high-resolution quasar spectra from the 
ESQ UVES/VLT archives (Dessaugcs-Zavadsky et al. 2003, 
hereafter Paper I). Sub-DLAs can be easily picked up since 
at these column densities damping wings are already formed 
and the rest equivalent width is Wrest > 2.5 A. 

In this paper, we analyse the sample of sub-DLAs pre- 
sented in Paper I. We first use these data to establish the 
column density distribution in a range unprobed before. 
These results have direct implications for the determination 
of the neutral gas content of quasar absorbers and its evolu- 
tion with redshift. We also emphasize the limitations of the 
present sample in terms of redshift path surveyed. The clus- 
tering properties of our sample of sub-DLAs are described in 
section 3. In section 4, we discuss the kinematic properties of 
the low- and high-ionization transitions associated with the 
sub-DLAs presented in Paper I. These results can be used to 
directly constrain some of the semi-analytic models recently 
presented in the literature. In section 5, we present the re- 
sults of chemical abundance determinations in sub-DLAs in 
conjunction with the abundances in DLAs taken from the 
literature and provide information on the global metallic- 
ity evolution of both individual quasar absorbers and the 
Universe. We also cross-correlate these properties with the 



kinematic information wo present in previous sections. Rela- 
tive abundances are then analysed and compared with recent 
models of galaxy formation and evolution, putting further 
constraints on the morphology of quasar absorbers. 

2 STATISTICAL PROPERTIES OF SUB-DLAS 

2.1 Sample Definition 

The sample of quasars used for this study is presented in 
Paper I and their characteristics are summarised in Table 1. 
This lead to a set of 22 quasars studied in various observa- 
tional programmes. In order to build the "statistical sam- 
ple" , we have ignored the quasars which have been targeted 
for the study of the Lya forest (for the analysis of voids, 
low column density column density distribution, etc). In- 
deed these have been preselected for not having a LLS in 
their spectrum, precisely the type of feature we are looking 
for in the present study. As reported in Paper I, a sub-DLA 
at Zabs ~ 1.838 was found in QllOl— 264, a quasar observed 
during UVES Science Verification. This absorber is included 
in the abundance analysis, although for consistency, neither 
the redshift path of the quasar nor the sub-DLA are used in 
the "statistical sample". This selection process resulted in 
a sample of 21 quasars suitable for studying the statistical 
properties of sub-DLAs. 

Similarly, sub-DLAs which do not have their Lya fear 
ture within the spectral coverage are automatically ex- 
cluded, although the determination of HI column density 
might be possible from other lines of the Lyman series. 
This is the case of the sub-DLA at Zahs = 3.171 towards 
Q1451-I-123, for which we undertook a detailed abundance 
study similar to the other systems, but which is not included 
in our "statistical sample". The resulting sample of sub- 
DLAs is thus composed of 10 systems. Although we looked 
for any system with N(HI) > 10^® atoms cm~^, the smallest 
column density sub-DLA detected in this sample is N(HI) 
= lO^^-^^ atoms cm~2. 

The large majority of the quasars in the sample were 
observed with UVES because there was a DLA along their 
line-of-sight, as can be attested from Table 1. There is no 
a priori reason that this creates a bias for our statistical 
analysis. 

Table 1 lists the minimum (zmin) and maximum (zmax) 
redshifts along which a sub-DLA could be detected along 
each quasar line-of-sight. Zmin corresponds to the point 
where the signal-to-noise ratio was too low to find absorp- 
tion features at the sub-DLA threshold of Wrest = 2.5 A, 
and Zma-x is 3000 km s~^ blueward of the Lya emission of 
the quasar. We took care to exclude the DLA regions and 
the gaps in the spectrum due to non-overlapping settings 
when computing the redshift path surveyed. 

2.2 Number Density 

The number density of quasar absorbers is the number of 
absorbers, n, per unit redshift dz, i.e., dn/dz = n{z). dz is 
computed by summing up the redshift paths surveyed along 
the line-of-sight of each of the quasar studied as given in 
Table 1. n{z) is a directly observable quantity, although, its 
interpretation is dependent on the geometry of the Universe. 
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Table 1. Az covered by each quasar used from the UVES archives and the associated sub-DLA systems. 
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" Zmin corresponds to the point where tiie signal-to-noise ratio was too poor to find absorption features 
and Zmax is 3000 km s~^ bluewards of the Ly-a emission of the quasar. The holes in the redshift intervals 
correspond to either a DLA or a gap in the spectrum due to non-overlapping settings. 

This quasar has been observed as part as the Science Verification of UVES for the study of the Lya 
forest. As such, it does not fulfill our criteria for its redshift path to be included in our "statistical sample" . 
Nevertheless, an analysis of the spectrum reveals the presence of a sub-DLA system at z^bs = 1.838 and we 

undertook the determination of the abundance of this absorber. 
This quasar is gravitationally lensed. Only the brightest line of sight was included in our study. 

This system is situated outside our spectral coverage. 

^ Our spectrum does not cover the Lya line of this sub-DLA (zabs = 3.171, Aobs ~ 5070 A). It has been 
identified thanks to the HI column density measurement reported by Bechtold (1994) and Petitjean, Sri- 
anand & Ledoux (2000). Several lines from the Lyman series and metals are available, hence, we undertook 
the determination of the metal content of this system but did not include it in our "statistical sample" . 



Indeed, the evolution of the number density of absorbers 

with redshift is the intrinsic evolution of the true number of 
absorbers combined with effects due to the expansion of the 
Universe. 

Measuring the incidence of sub-DLAs down to N(HI) 
= 10^^ atoms cm~^ in a given sample of quasars requires 
a number of high resolution spectra in order to unambigu- 
ously select all the absorption systems. Nevertheless, indi- 
rect information is also provided by the number of LLS in 
a given quasar sample. Peroux et al. (2003) fitted the ob- 
served cumulative number of absorbers of both DLAs and 



LLSs with a F-distribution as suggested by Pei & Fall (1995) 

and first implemented by Storric-Lombardi, Irwin & McMa- 
hon (1996a). This puts constraints on the number of systems 
in the 10^'' < N(HI) < 2 x 10^° atoms cm"^ column density 
range. They then used these predictions to compute the ex- 
pected number density of sub-DLAs. The results from this 
work are tabulated in the last column of Table 2 for com- 
parison with the observations presented here. 

Excluding QllOl— 264, we have sampled a total redshift 
path dz = 17.5 and find a total of 10 sub-DLAs. The com- 
parison with the predictions of Peroux et al. (2003) show 
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Table 2. This table gives the redshift path surveyed, dz, and 
corresponding physical distance interval, dX (in a Q,m = 0.3, 
Qa = 0-7 cosmology) for various redshift ranges. The observed 
number density of sub-DLAs is from our "statistical sample" of 
21 quasars from the UVES archive spectra. The predicted number 
density is from Peroux et al. (2003) who used a fit to the observed 
number of DLAs with an additional constraint from the total 
number of LLS. 



z range 


dz 


dX 


number 


n{z) 


n{z) 




obs 


obs 


of sub-DLAs 


obs 


predicted 


<3.5 


14.4 


47.0 
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0.49 


0.46 


>3.5 
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1.66 


1.5-2.5 
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1 


0.53 


0.67 


3.5-4.0 


1.7 


6.8 


2 


1.18 


1.52 


4.0-5.0 


1.4 


5.7 


1 


0.71 


2.00 


total 


17.5 


59.5 


10 


0.57 


0.58 



a good agreement in the intermediate redshift range. At 
z > 3.5, the low number statistics makes the results more 
uncertain. Nevertheless, it is important to note that in the 
quasar sample studied here all but 4 of the objects have 
Zcm < 3.5. The sub-DLAs are thus predominantly located 
in the intermediate redshift range 2 < Zabs < 3.5. Indeed, 
we find 7 sub-DLAs at Zabs < 3.5, where we have surveyed 
dz — 14.4, while we find 3 sub-DLAs at Zabs ^ 3.5, where we 
have only surveyed dz = 3.1. Consequently, the determina- 
tion of n{z) at Zabs 3.5 is not possible with the current set 
of data. The number of sub-DLAs and associated redshift 
paths are summarised in Table 2 for various redshift ranges. 

2.3 Column Density Distribution 

2.3.1 Observationally Constraining f{N) down to log 
N{HI) = 19.0 

Our sample of sub-DLAs from the UVES archive quasar 
spectra can be further used to determine the column den- 
sity distribution of absorbers in the column density range 
10" < N(HI) < 2 X 10^° atoms cm"^ The column density 
distribution describes the evolution of quasar absorbers as 
a function of column density. It is defined as: 



f{N, z)dNdX = 



■dNdX 



(1) 



where n is the number of quasar absorbers observed in 
a column density bin [A*', A*' + AN] obtained from the obser- 
vation of m quasar spectra with total absorption distance 
coverage X/I^i ^-^i- The distance interval, dX, is used to 
correct to co-moving coordinates and thus depends on the 
geometry of the Universe. For a non-zero A-Universe: 



Xiz) 



(l+zf [(1 + zf{l + zQm) - z(2 + z)nA] 



-1/2 



dz{2) 



The values of dX used in our study are given in Table 
2 for an Qm = 0.3, Qa ~ 0.7 cosmology. We have sampled a 
total redshift path dz = 17.5. This corresponds to a distance 
interval dX = 59.5 in Qm = 0.3, Qa = 0.7 cosmological 
model. 
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Figure 1. The solid bins at log N{HI) > 20.3 atoms cm are 
measurements of the column density distribution of DLAs from 
Peroux et al. (2002). The dashed bins (17.2 < log N{HI) < 20.3 
atoms cm~^) are deduced from the fit to the observed cumula- 
tive number of quasar absorbers (Peroux ct al. 2003), while the 
solid bin in the range 19.0 < log N{HI) < 20.3 atoms cm~^ 
is the observed number density of sub-DLAs from the UVES 
archive quasar spectra presented in this paper. Models (corrected 
to flui = 0.3, Qa = 0.7 cosmology) of Zheng & Miralda-Escude 
(2002) with halo masses of lO^^, lO", 10^0, 10^ Mq (from top to 
bottom) are overplotted. The vertical normalisation of the curves 
is unconstrained in the models and thus on this figure the lines are 
arbitrarily shifted to best fit the data points. On the other hand, 
the shape of the model distributions can be directly compared 
with the data in order to determine the neutral fraction in the 
absorber at the radius where self-shielding starts to be important 
(see text for further explanation). 



We have used the sample of sub-DLAs from the UVES 
archive quasar spectra to determine the column density dis- 
tribution down to N(HI) = 10^^ atoms cm~^ for the first 
time. The results are shown in Figure 1 together with the 
predictions of Peroux et al. (2002) computed from the ex- 
pected number of sub-DLAs and a compilation of DLAs at 
all redshifts. The observational results presented here are 
consistent with the predicted change of slope of the col- 
umn density distribution at N(HI) ~ 10^^ atoms cm~^, as 
suggested by previous work (Petitjean et al. 1993; Storrie- 
Lombardi & Wolfe 2000; Peroux et al. 2003). This feature is 
probably the signature of the onset of self- shielding in these 
systems. 



2.3.2 Comparison with Models 

Although simulating high column density systems such as 
DLAs is still extremely challenging because of resolution 
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limitations, Gardner ct al. (1997) have tried to overcome 
the problem by imposing the density profile of resolved halos 
onto unresolved ones. They derive the evolution with red- 
shift of the column density distribution and, interestingly, 
they find a flattening of the distribution somewhere in the 
region log N (HI) = 18.5 — 20.0 atoms cm~^. Nevertheless, 
the theory predicts little change in the form of the column 
density distribution function over the redshift range 2 — 4. 
This seems counter to current observations since a strong 
evolution of the rmmbor density of quasar absorbers with 
redshift has been derived (Peroux et al. 2002), implying a 
factor of 3 or more difference in the number of LLS between 
z = 2 and z = A. 

More importantly, the new set of observations can be 
directly compared with models of the column density distri- 
bution (Corbelli, Salpeter & Bandiera, 2001 and Zheng & 
Miralda-Escude, 2002). Zheng & Miralda-Escude (2002), in 
particular, modeled DLAs as spherical isothermal gaseous 
halos ionized by external cosmic background to predict the 
column density distribution. Their models for different halo 
masses: 10^^, lO", 10^°, lO'' Mq (from top to bottom) are 
overplotted in Figure 1. The normalisation of the curves is 
arbitrary since it depends, among other things, on the size 
distribution of clouds, which is not known, but the shape can 
be directly compared with the data. In addition, Zheng & 
Miralda-Escude show that the column density at which the 
flattening takes place depends only on the neutral fraction in 
the absorber at the radius where self-shielding occurs. The 
present data exclude the models with high mass halos. 

In addition, according to Peroux et al. (2002), the 
change of slope of the column density distribution should 
evolve with redshift, moving out in redshift from z = 2 to 
2 = 4. Unfortunately, the current sub-DLA sample is not 
sufficient to study the evolution with redshift of the distri- 
bution down to N(HI) = 10^® atoms cm~^. More data are 
required at high redshift in order to be able to probe the 
evolution with time of the properties of the sub-DLAs. 



3 CLUSTERING OF SUB-DLAS 

As already noted all the quasars studied here were extracted 
from the ESO UVES archive and so the large majority were 
observed because there was a DLA along their line-of-sight. 
Although the DLAs themselves were excluded from the red- 
shift path surveyed in the quasar, it is of interest to note 
that on at least one occasion, the sub-DLA discovered is 
situated in the wing of a known DLA: 2sub-DLA = 2.557 
and 2DLA = 2.466 towards Q1223+1753, giving Az = 0.09, 
corresponding to Aw ~ 8000 km s^^ (see Paper I). In addi- 
tion, we note that out of the 12 sub-DLAs which make up 
our sample, at least 3 systems (^^sub-DLA = 1.996 towards 
Q2116-358, «sub-DLA = 3.565 towards PSS J2155-M358 
and Zsub-DLA = 3.882 towards PSS J2344-I-0342) cannot 
be fitted with only one absorption system, suggesting the 
presence of smaller column density absorbers in their close 
vicinity (see Paper I for more details on these systems). Ob- 
viously, these features illustrate the difficulty in determin- 
ing the HI column density of sub-DLA systems and differ- 
entiating these objects from a conglomerate of smaller col- 
umn density objects. In the particular cases of ^sub-DLA = 
3.356, 3.589 and 4.214 towards BR J0307-4945, we cannot 



estimate the size and number of absorbers involved. We refer 
the reader to section 3.1 of Paper I for discussion of these 
specific cases. 

To summarise, our sample of 12 sub-DLAs contain 4 
systems associated with other absorbers, one of which is 
a DLA, the remaining being smaller column density sys- 
tems. Erom Peroux et al. (2002), the expected probability 
of absorbers with N(H) > 10^^ ° atoms cm"^ is n{z) ~ 0.4 
at 2 ~ 2.5, if the absorbers where randomly distributed. 
Thus there is only a 1 in 50 chance of finding quasar ab- 
sorbers separated by Az = 0.09 {/^v 8000 km s"^). The 
observed incidence of clustered absorbers (1 out of 12) is 
therefore higher than expected. For comparison, out of the 
hundreds of DLAs known today, only 3 cases of such clus- 
tering along a given line-of-sight have been reported in the 
literature. Prochaska & Wolfe (1999) report a pair of DLAs 
towards Q2359— 02 at 2 = 2.1 and separated by At; ~ 5700 
km s-i. Ellison & Lopez (2001) have jirescnted a pair of 
(strictly speaking) DLA/sub-DLA at z = 1.8 separated by 
Aw ^ 2000 km s~^ which are characterised by a low o/Fe 
elemental ratio. Lopez et al. (2001) have found 3 DLAs at 
z = 2.6 in CTQ 247 separated by Av ~ 6000 km s"^ which 
show similar abundance patterns (Lopez & Ellison 2003). 

The quantitative interpretation of this phenomenon is 
not straightforward. Lyman Break selected galaxies (e.g. 
Steidel et al. 1998) and Ly-a emitting galaxies (Ellison et al. 
2001) have been shown to be clustered at high redshifts. Sim- 
ilarly, studies of transverse clustering from quasars pairs or 
groups show a highly significant overdensity of strong ab- 
sorption systems over separation lengths from ~ 1 to 8 h~^ 
Mpc (D'Odorico, Petitjean & Cristiani 2002). Our findings 
are therefore not surprising in that respect. In all cases, pairs 
of sub-DLAs/DLAs are potentially of great interest since, in 
such cases, one could reasonably assume that the incident 
radiation field is similar for both absorbers, thus provid- 
ing further information on the otherwise poorly constrained 
photoionization spectrum. Such studies would nevertheless 
require a good estimate of HI column density of each sys- 
tem (often requiring other lines down the Lyman series) and 
that the metals associated with each absorber are well char- 
acterised. 



4 KINEMATICS OF SUB-DLAS 

The kinematic properties of quasar absorbers can be used 
to provide further information on their nature. Using var- 
ious line diagnostics to parameterise the symmetry, veloc- 
ity width and the edge-leading profiles of the metal lines of 
17 DLA systems, Prochaska & Wolfe (1998) suggest that 
the absorbers are fully formed, large, rapidly rotating galac- 
tic discs with Vcirc ^200 km s~^. More recently, Wolfe & 
Prochaska (2000b) have used a sample of 35 DLAs (Wolfe & 
Prochaska 2000a) and compared them with semi-analytical 
models of galaxy formation. They note that none of the mod- 
els reproduce the overlap between low- and liigli-ionization 
metal profiles observed in the data. In addition, this scenario 
is in stark contreist with the currently favoured hierarchi- 
cal structure formation models, where present-day galaxies 
are assembled from virialised sub-units over a large redshift 
range (2 ~ 1 — 5). Other hydrodynamic N-body simulations 
and semi-analytical models (Haehnelt, Steinmetz & Rauch 
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Figure 2. Velocity space superposition of tlie low- (Fe II or Si ll) and liigli-ionization (C IV or Si iv) transition profiles for the sub-DLAs 
of our sample where high-ionization transitions are detected. The two states of the gas overlap in most cases apart from x^bs = 3.171 in 
Q145H-123. 
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Figure 3. Comparison of the kinematic properties of DLAs (black line) and sub-DLAs issued from both our sample and the literature 
(light coloured histogram). The parameters plotted are (from left to right, from top to bottom) Aio„, Ajjigh, 6v and /ratio- Although 
the statistics on sub-DLAs are still small it seems that their properties roughly span the parameter space of the DLAs, thus challenging 
multicomponcnt semi-analytical models. 



1998; Mailer et al. 2001; Mailer at al. 2002) have shown that 
the low ion absorption profiles can equally well be inter- 
preted as a signature from merging proto-galactic clumps in 
collapsing dark matter halos with small virial velocities (~ 
100 km s~^). In addition, Ledoux et al. (1998) have used a 
sample of 26 DLAs and find a correlation between the asym- 
metry and AV for AV 150km s"'^. They suggest that this 
correlation is evidence that rotation velocities may dominate 
the narrower metal lines. 



4.1 Sub-DLAs Characteristics 

First we examine the number of components necessary to 
fit the low-ionization transitions of the sub-DLAs, noting 
that out of 12 systems, 5 require more than 10 components. 
We also note that most of the remaining systems (another 
5 absorbers) are well fitted with only 3 or less components. 

We chose to describe the properties of the lines of the 
low- and high-ionization transitions in our sample of sub- 
DLAs with the help of four simple quantities, namely: 



(i) Alow: the width of the low-ionization transitions as 
measured from the most saturated Fe ll or Si ll lines. 

(ii) Ahigh: the width of the high-ionization transitions as 
measured from the most saturated C iv or Si iv lines. 

(iii) 5v: the offset between the mean velocity of the high- 
and low-ionization transitions. 

(iv) /ratio: the ratio of the velocity widths of the high- 
and low-ionization transitions. 

All these parameters are used to quantify the differ- 
ences between various states of the gas observed in absorp- 
tion. The shift between the low- and high-ionization tran- 
sitions for our sample of sub-DLAs is shown in Figure 2 
for the 8 systems where high-ionization transitions are de- 
tected. In most cases, although the two sets of absorption 
overlap, the profiles are widely different. In one particular 
object (zabs = 3.171 in Q1451-I-123), the Si ll is not super- 
imposed at all on the C iv doublet in velocity space, the 
two being separated by ~ 75 km s~^. It might be that the 
C IV doublet is actually not associated with the sub-DLA 
itself (see Paper I for more details). As in DLAs, the high- 
ionization transitions spread over larger velocity intervals 
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than their low-ionization counterparts, but their correspon- 
dence suggests that the systems are within the same poten- 
tial well (Wolfe & Prochaska 2000b). However, in two cases 
(«abs = 3.882 in Q2344+034 & z^bs = 2.087 in Q1444+014), 
the low-ionization transitions traced by Si ll spread over a 
larger velocity interval than the high-ionization transitions. 



4.2 Comparison with DLAs 

In order to compare the kinematic properties of sub-DLAs 
with DLAs, we have determined the parameters described in 
the previous section for 72 DLAs taken from the literature. 
The systems quoted are from a mixture of UVES/HIRES 
and lower resolution spectrographs for high-rcdshift targets, 
and from HST at lower redshifts. Thirteen additional sub- 
DLAs are found in the literature. We made the relevant 
measurements from the quasar spectra whenever we had 
them, otherwise we directly used the original publications 
to quantify the width of the ion lines. We preferentially used 
the Fc II and C iv linos whenever possible. In addition, if the 
data came from the UCSD HIRES /Keckl survey (Prochaska 
& Wolfe 1998, 2000 and 2002; Prochaska, Gawiser & Wolfe 
2001; Prochaska et al. 2002), we took the velocities used for 
the Optical Depth Method given on the survey's web page ^ 
kindly provided to the general community and maintained 
by Jason X. Prochaska. 

The resulting distributions for Aiow, Ahigh, Sv and /ratio 
are illustrated in Figure 3. Although the number of sub- 
DLAs is still small, it appears that their velocity spreads 
in both low- and high-ionization transitions are statistically 
equivalent to the DLAs (top left and top right panel of Fig- 
ure 3). In most cases, the shift between the low- and high- 
ionization transitions is small in DLAs (less than 20 km s~^) 
(bottom left panel of Figure 3) . Finally, the /ratio values of 
DLAs and sub-DLAs are comparable (bottom right panel 
of Figure 3). These measurements are important to con- 
strain semi-analytical models. Indeed, Mailer et al. (2002) 
find that a large fraction of absorbers with log N(HI)> 20.0 
cm~^ might be composed of multiple component discs. In 
contrast, they predict that the smaller systems (in particu- 
lar those with 4 x lO"^ < N(HI)< 10^" cm"^) will mostly be 
single discs. Therefore, Sv is expected to be smaller in these 
systems than in classical DLAs. Although our data set is still 
small, the observations presented here do not support this 
interpretation, thus challenging multiple component models. 
A possible explanation is that some of the low ionization gas 
in sub-DLAs is in hot gas. Indeed, the data show a hint of 
decreasing fratio with lower column densities, therefore sug- 
gesting that the high and low ions are in the same gas at 
lower N(HI). 

5 CHEMICAL ABUNDANCES OF SUB-DLAS 
5.1 Ionization Fraction 

In Paper I, we have presented the hydrogen and ionic column 
densities of all of the sub-DLAs included in Table 1, using 

■t http: / /kingpin.ucsd.edu/ ~hiresdla/ 



Table 3. Results from the order ranking Kendall test and slopes 
of tiie fits assuming a linear correlation for different sub-sets of 
quasar absorbers, ignoring all upper and lower limits. The lack of 
correlation between the redshift and [Fe/H] is due to the larger 
column density systems (log N(HI) > 21.0). In contrast the sub- 
DLAs appear to show more correlation, with a more pronounced 
slope, although the results are still weakened by low number 
statistics. 



data set nber of t c.l. Slope 

definition abs 

log N(HI) > 20.3 72 -0.29 99.53% -0.18 ± 0.12 

log N(H1) < 20.3 17 -0.50 84.27% -0.40 ± 0.22 

log N{HI) < 21.0 75 -0.40 99.99% -0.28 ±0.11 

All AlisorlxTs 89 n.:-!5 !)<).r,:-!'X -(1,21 ±0.11 



standard Voigt profile fitting. Since their relatively low col- 
umn density implies that some of the gas might be ionized, 
the ionization state of the sub-DLAs has been investigated 
in detail in Paper I using the photoionization models of the 
CLOUDY software package. The determination of the ion- 
ization parameter U for a given absorber system crucially 
depends on the detection of intermediate-ionization transi- 
tions. In addition, the photoionization model relies on sev- 
eral input parameters related to the nature of the absorber 
and the characteristics of the radiation field. In this context 
we chose to extract from our analysis a global ionization 
correction rather than individual corrections limited by par- 
tial information. We found the correction to be < 0.3 dex 
for most elements down to N(HI)= 19.3 cm^'^. We also note 
that the correction to apply is quite different from species to 
another. In the case of O l or C II, for example, the ioniza- 
tion fraction is small, while it becomes important for such 
elements as Al ii or Zn ii. One needs to keep these facts in 
mind while interpreting the various figures presented in the 
following sections. 

5.2 Metallicity Evolution 

In this section, we present the results of the metallicity study 
of the sample of sub-DLAs and compare them with the char- 
acteristics of the well studied DLAs. The references for ab- 
sorbers from the literature are the same as in Section 4.2. 
The chemical abundances in DLAs have been used to trace 
the metallicity evolution of gala^Kies over a large look-back 
time (e.g. Pettini et al. 1999; Prochaska & Wolfe 2000) and 
contrary to virtually all chemical models, the observations 
in DLAs indicate mild evolution with redshift. By extend- 
ing the abundance determination to systems of lower column 
densities, such as the sub-DLAs, we can gain insight in the 
evolution of the galactic chemical evolution of galaxies, since 
they are more likely to represent the basic building blocks 
of hierarchical growth of structure. 

5.2.1 Metallicity of Individual Systems 

[Fe/H]: Given that the ionization correction on the Fe ii 
column density is within the observational errors, the abun- 
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Figure 4. [Fo-Zn/H] mctallicitios for DLAs (squares) and sub-DLAs (open circles for data from the literature, filled circles for the 
sub-DLAs from this work) as a function of redshift. The symbol sizes are proportional to the HI column density of the systems. The 
evolution with redshift of the [Fe/H] ratio is more pronounced for sub-DLAs than for DLAs (solid and dashed lines, respectively). The 
dashed-dotted (dotted) curves are results from models from Hou, Boissier & Prantzos (2001) for sub-DLAs (DLAs). 
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dance measurements of the sub-DLAs can be directly com- 
pared with those of the DLAs gathered from various sources 
in the hterature. We use the abundance ratios with re- 
spect to solar values defined in the usual manner § and 
presented in Paper I. For the sake of homogeneity wc re- 
compute these ratios for all data from the literature assum- 
ing [Fe/H]Q = -4.50 (Grevesse & Sauval 1998). The top 
panel of Figure 4 presents the evolution with redshift of 
the [Fe/H] absolute abundances in individual systems. The 
symbol sizes in this figure arc proportional to the HI column 
density of the systems. The classical DLAs are represented 
with square symbols while the sub-DLAs are circles (open 
for data from the literature, filled for the data from the sam- 
ple presented in this paper). Wc applied the order ranking 
Kendall test to several subsets of the sample as shown in 
Table 3, ignoring all upper and lower limits. As previously 
reported (e.g. Prochaska & Wolfe 2002; Dessauges-Zavadsky 
et al. 2001), there is evidence for a mild evolution of the in- 
dividual [Fe/H] in DLAs (correlation coefficient of -0.29 at 
99.53% confidence level). The analysis suggests a stronger 
correlation of the [Fe/H] in sub-DLAs although the signifi- 
cance of the test is weakened by low number statistics (cor- 
relation coefficient of 0.50 at 84.27% confidence level for a 
total of 17 sub-DLAs). The effect is predominant ai z < 2 
and most of the evolution observed lies in this redshift range. 
We fitted the two populations of absorbers assuming a linear 
correlation and found a slope of a = —0.18 ± 0.12 for DLAs 
and Q = —0.40 ±0.22 for sub-DLAs (dashed and solid line of 
the top panel of Figure 4 respectively). The evolution with 
redshift of the [Fe/H] ratio might be more pronounced for 
sub-DLAs than for DLAs. This difference, if confirmed by a 
larger sample of data, suggests that the detection of DLAs 
is more biased by dust at low redshift, or that the sub-DLAs 
are associated with a class of gala^Kies which better trace the 
overall chemical evolution of the Universe. 

[Zn/H]: The well-known drawback of using Fe ll for 
abundance determination is the fact that this element is 
sensitive to depletion onto dust grains (Pettini et al. 1997). 
An alternative to overcome the dust depletion issue is pro- 
vided by Zn, an element known to be only mildly depleted. 
The lower panel of Figure 4 present the evolution of the 
[Zn/H] with redshift. Again for the sake of homogeneity we 
recompute these ratios for all data from the literature as- 
suming [Zn/H]Q = —7.33 (Grevesse & Sauval 1998). In our 
sample of sub-DLAs, most of the Zn ii lines are weak or un- 
detected and thus wc can only provide upper limits. In addi- 
tion, the ionization correction in sub-DLAs is important for 
Zn 11 and thus the interpretation of the abundances is more 
complicated. We again used the order ranking Kendall test 
for the subsets of the Zn ii sample. When considering both 
DLAs and sub-DLAs (and excluding all limits), we deter- 
mine aslope a = — 0.28±0.23, i.e. slightly shallower than the 
one derived by Vladilo et al. (2000) and Dessauges-Zavadsky 
et al. (2001) (a ~ —0.3) since the latter works loosely in- 
clude sub-DLAs in their samples. Nevertheless, the signif- 
icance of this result is low (correlation coefficient of -0.25 
at 84.27% confidence level for a total of 40 abundances). 
Analysing DLAs only, leads to a slope a = —0.21 ± 0.25 



§ [X/H]=log[N{X)/N{H)]DLA-iog[N{X)/N(H)]Q, assuming 
that N{X) = N{XII) and N{H) = N{HI). 



for 36 systems. There are not enough Zn ll measurements 
in the sub-DLAs range to allow such analysis. An impor- 
tant point, however, is the fact that the slope in the Zn ii 
DLA sub-set (a = —0.21 ±0.25), a tracer of true metallicity 
evolution, is not as steep as in the Fc n sub-DLA sub-set 
(a = —0.40 ± 0.22). This suggests that the evolution we 
probe in these latter absorbers is not due to the hidden ef- 
fect of dust. 

Comparison with Models: These observational re- 
sults can be compared with models of galaxy evolution and 
the corresponding evolution of metal abundances to fur- 
ther test the nature of quasar absorbers. Complications arise 
from the fact that the observed systems might have different 
formation epochs, a parameter which cannot be directly con- 
strained by the data. Prantzos & Boissier (2000) use models 
of the chemical evolution of disc galaxies, calibrated on the 
Milky Way (Boissier & Prantzos 1999, 2000), which they 
then compare with the Zn abundances of DLAs. This com- 
parison is extended to other elements in Hou, Boissier and 
Prantzos (2001). They find that the results are compati- 
ble with available observations, provided that several bi- 
ases are taken into account, one of these biases being re- 
produced by applying a "dust filter" to the models. This 
filter was proposed by Boisse et al. (1998) who suggest that 
high column density and high-metallicity systems are not 
detected because the light of background quasars is severely 
extinguished. This claim arises from the observed trend in 
[X/H] versus log N(HI) in quasar absorbers (sec section 5.2.3 
for more on this issue). The results of the models are fur- 
ther constrained by the HI column density detection lim- 
its in current surveys. In Hou, Boissier & Prantzos (2001), 
this limit corresponds to the formal definition of DLAs, 
logN(HI) > 20.3. The mean [Fe/H] value of the model is 
presented in Figure 4 as a dotted curve. The scatter in the 
observations is genuine (since it is higher than the error bars 
on individual metallicity measurements) and is well repro- 
duced by the model (sec Figure 2 of Hou, Boissier & Prant- 
zos 2001). 

Adjusting the threshold to the sub-DLA definition re- 
sults in the dashed-dotted curve presented in Figure 4. The 
models do not indicate a change in evolution with redshift 
between DLA and sub-DLA systems (see dotted and dashed- 
dotted curves in the bottom panel of Figure 4) . Rather they 
suggest that the metallicity in sub-DLAs is lower than in 
DLAs. This is a natural consequence of the abundance gradi- 
ents observed in nearby galaxies and in the Milky Way, well 
reproduced by the models (Prantzos & Boissier 2000). On 
the contrary, the observations indicate that the sub-DLAs 
seem more metal-rich than the DLAs, at least at low red- 
shift. Therefore although the disc models are in quite good 
accordance with the DLAs observations, the mismatch be- 
tween these models and the sub-DLA observations suggests 
that sub-DLAs do not arise in disc galaxies. 

The [Zn/H] observed redshift evolution can also be com- 
pared with models of galaxy evolution proposed by Hou, 
Boissier & Prantzos (2001). It should be emphasized that 
in the models, the Zn ii abundance is set to follow the 
Fe abundance. The difference between the two arises only 
from different "biases"; the current Zn detection can be 
accounted for by applying to the models the filter: log 
N(HI)-|-[Zn/H]>18.8, an estimate of the Zn limit detection 
(e.g. Prantzos & Boissier 2000). Again, setting the column 
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Figure 5. HI column density-weighted mean metallicities for various sub-sets of quasar absorbers. The dotted bins are for constant HI 
intervals and the solid bins are for constant redshift intervals. The evolution of [(Fc/Hdla)] is clearly more pronounced for sub-DLAs 



than for DLAs, but this result is not apparent when all absorbers with N(HI) > 10^^" 
the number of DLAs (72) is much larger than the number of sub-DLAs (17). 



atoms cm are considered. Indeed, in our sample 
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density cut-off in the models to the DLA definition repro- 
duces the trend well (dotted curve in the bottom panel of 
Figure 4) and the scatter between measurements (see also 
Prantzos & Boissier 2000). Once the cut-off is chosen to 
match the sub-DLA definition (dasliod-dottcd curve in bot- 
tom panel of Figure 4), the models suggest that the Zn ii 
abundances in the systems at redshift greater than ~ 1.3 
are below current detection limits. This is indeed what is 
observed since most of the [Zn/H] presented here are upper 
limits. Note that to be able to measure the Zn abundance 
(log N(HI)-|-[Zn/H]>18.8) in sub-DLAs, the system must 
satisfy : [Zn/H] > -1.5. 

An alternative modelling approach has been recently 
proposed by Ceri et al. (2003) using hydrodynamic simu- 
lations to reproduce the gas and metal evolution observed 
in DLAs. They find that the slow metallicity evolution cur- 
rently observed in DLAs can be explained by the sequential 
formation of galaxies from the highest metallicity systems. 
This effect would be combined with the increase of metallic- 
ity of the absorbers with the smallest abundances. Including 
an artificial cut-off in order to mimic the presence of dust 
in DLAs, the simulations are found to be in excellent agree- 
ment with the observations (Gen et al. 2003). 



5.2.2 Weighted Mean Metallicity 

Since DLAs and sub-DLAs are known to contain a major 
fraction of the neutral gas at all redshifts (Peroux et al. 
2003), all of those systems should be used to probe the 
metallicity of the Universe and its evolution with time. A 
quantitative way of estimating this evolution is provided by 
the HLweighted mean metallicity. Pettini et al. (1997) state 
that "under the working assumption that DLAs account for 
most of the material available for star formation at high red- 
shift, the quantity [(X/Hdla)] (X=Fe or Zn) is a measure 
of the degree of metal enrichment by the Universe at a given 
epoch" . 

Figure 5 presents the column density-weighted mean 
Fe 11 abundances of n systems in redshift bins of Az = 1 
(solid bins) and in constant HI interval (dashed bins), for 
various sub-samples of quasar absorbers. Although more 
advanced statistical methods have been recently suggested 
(Kulkarni & Fall 2002), we choose to follow the prescription 
from Pettini et al. (1997): 

[(Fe/Hj^i^^)] = log((Fe/H)j3i^^) - log (Fe/H)^ (3) 
where 

n 

((Fe/H)^^^) = ^ (4) 

E ^(Htotal)i 

and the error are estimated from the standard devia- 
tion, a': 

= (rA^'^l^^i - [<Fe/HoLA)l)'^ /(" - 1) (5) 

with 

[Fe/H.] = log(Fe/H). - log (Fe/H)^ (6) 



The results of these calculations are shown in Figure 5 
for different subsets of quasar absorbers. The top panel in- 
cludes all quasar absorbers (DLAs -|- sub-DLAs) as well as 
systems with N(HI) > 10^^ ° atoms cm~^. The bottom panel 
shows DLA and sub-DLAs separately. We note that the evo- 
lution of [(Fe/HDLA)] is clearly more pronounced for sub- 
DLAs than for DLAs, but this result is not apparent when 
all absorbers with N(HI) > 10^® ° atoms cm~^ are consid- 
ered. Indeed, in our sample the number of DLAs (72) is much 
larger than the number of sub-DLAs (17). On the contrary, 
in the Universe we expect the number of sub-DLAs to be 
up to 3 times (at x ~ 4) the number of DLAs (Peroux et al. 
2002). Since high column density systems dominate the HI- 
weighted mean metallicity, it will be important to increase 
the sub-DLA sample size to better probe [(Fe/HDLA)]. Nev- 
ertheless, these results show that the HI column density- 
weighted mean metallicity Fe ii of sub-DLAs do evolve with 
redshift more markedly than for the DLA population. The 
dotted bins correspond to constant HI intervals and there- 
fore present similar number of systems. They also reveal an 
increase of metallicity with decreasing redshift, showing that 
this result is not an artifact of the low number of sub-DLAs 
known at 2; < 2. Starting from a metallicity 1/100 solar at 
z ~ 4.5, the sub-DLAs evolve up to a metallicity 1/3 so- 
lar at z ~ 0.5. The ionization correction cannot explain the 
evolution observed since we have shown in Paper I that for 
[Fe ii/HI], it does not exceed 0.2 dex. These results again re- 
inforce the hypothesis that sub-DLAs better trace the global 
metallicity evolution and therefore should be included in the 
metallicity determination to obtain a complete picture of the 
abundance. 

5.2.3 Correlating with Other Properties 

[Fe-Zn/H] versus HI: Figure 6 presents [Fe-Zn/H] as a 
function of HI column density for DLAs and sub-DLAs in 
three different redshift ranges: z<2, 2<z<2> and « > 3. 
In the top panel ([Fe/H]), the spread of values is larger at 
lower HI column densities. This suggests that we are sam- 
pling a different class of objects, or at least different parts 
of the same object as suggested by e.g. Boissier, Peroux & 
Pettini (2002). Another possibility is that dust plays an im- 
portant role in the diversity of properties through a range 
of star formation histories. In contrast, there is little spread 
in the [Fe/H] ratio at very high HI column densities (log 
N(HI)> 21) posing the question of why these systems dif- 
fer from the other quasar absorbers. This might indicate 
that systems with high column densities belong to a well- 
defined population of objects, though more probably, the 
lack of apparent spread could be the result of a bias against 
high column density, high metallicity absorbers due to their 
high dust content (Boisse et al. 1998). Indeed, the bottom 
panel of Figure 6 shows that the observational bias noted by 
Boisse et al. in the HI versus [Zn/H] relation, still holds (i.e. 
no data are situated above the [Zn/H]-|-logN(HI) >21 solid 
line). A direct dust diagnostic in sub-DLAs is not straight- 
forward to make since Zn 11 is seldom detected. Using the 
three Zn II detections in our sample of 12 systems, we find 
that the [Zn/Fe] value is around 0.3-0.4 dex. This suggests 
that the amount of dust is not negligible, although naturally 
there are too few detections to draw any firm conclusions. In 
addition, a study of absorbers derived from a radio-selected 
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Figure 6. [Fe/H] and [Zn/H] as a function of HI column density in DLAs and sub-DLAs. The dotted line marks the DLA formal 
definition. On the bottom panel, the solid line corresponds to the "dust filter" proposed by Boisse et al. (1998), namely [Zn/H]+logN(HI) 
>21. 
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Figure 7. [Fe/H] as a function of various kinematic parameters: Aio„ and Aj^igh (see Section 3 for details). The symbol shapes depict 
different redshift ranges. The open symbols are for DLAs and the filled symbols are for sub-DLAs. The boxes represent the mean in a 
given velocity interval with rms errors and suggest an increase of metallicity towards larger ionization widths. 



quasar sample (CORALS) does not indicate that a dust bias 
greatly afltects DLAs studies (Ellison et al. 2001), but in the 
authors own words the conclusions are only tentative due to 
the limited size of the data set. Furthermore, in the range 
20.0 <log N(HI)< 21.0, i.e. the bulk of the data available, 
the high redshift objects are less evolved than the others. 
In any case, we find that the metallicity evolution of quasar 
absorbers is a strong function of the HI column density. In 
addition, we have shown that the ionization correction in 
Fe 11 (see Paper I) is similar within a group of systems with 
similar column densities. In the sub-DLA range, it tends to 
increase the Fe ll column density measured by 0.1-0.2, thus 
not affecting the apparent evolution. 

[Fe-Zn/H] versus kinematic properties: Figure 7 
presents the metallicity of quasar absorbers as a function of 
their kinematic properties (Aiow, Ahigh) as measured in Sec- 
tion 3. The symbol shapes denote different redshift ranges. 
The open symbols are for DLAs and the filled symbols are 
for sub-DLAs. The results indicate that in most cases the 
properties are homogeneously mixed between samples of dif- 
ferent HI column densities or redshifts. The boxes represent 
the mean in a given velocity interval with rms errors and 
suggest an increase of metallicity towards larger ionization 
widths, although the statistical significance of this result is 
low. 

5.3 Relative Abundances 

The comparison of elemental ratios, independent of HI col- 
umn density, but as a function of metallicity is also very 
important. We concentrate in particular on the a (i.e. O, 
Si) over Fe-peak element ratios versus [Fe/H] metallicities, 
since this relation provides information on the star forma- 
tion histories of quasar absorbers. [a/Fe] versus [Fe/H] is a 
strong function of the star formation history: a/Fe depends 
on the lifetimes of the element progenitors, whereas [Fe/H] 



depends on the star formation rates. The a-elements are 
typical products of Type II supernovae, which are likely to 
dominate the early stages of galaxy formation. Fe-peak ele- 
ments on the contrary are products of Type la supernovae, 
which play a role in later evolutionary stages. Therefore, 
high overabundances of a-elements at low [Fe/H] indicate 
recent star formation, while solar a-element abundances at 
low [Fe/H] suggests a lower star formation rate or bursts 
separated by quiescent periods. Nevertheless, the interpre- 
tation of these elemental ratios is once again complicated 
by differential depletion: Fe ll and Si ll are known to be de- 
pleted onto dust grains in different proportions while O I is 
expected to remain unaffected by dust depletion. 

[Si/Fe]: Figure 8 shows that [Si/Fe] is similar in both 
DLAs and sub-DLAs, averaging around -1-0.5 dex. For the 
sake of homogeneity we recompute these ratios for all data 
from the literature and models assuming [Si/H]Q = —4.44 
(Grevesse & Sauval 1998). This value is typical of Galactic 
metal-poor stars and is a strong indicator of a enhancement 
if dust depletion does not play a major role. Nevertheless, 
the differential depletion of these two elements can lead to 
the same enhancement as the one observed. This observation 
can be directly compared with various models of galaxy for- 
mation and evolution. Models of Hon, Boissier & Prantzos 

(2001) , with cut-off at the DLA and sub-DLA column den- 
sity definitions, are presented in Figure 8 as dark-coloured 
dashed and solid contours, respectively. In both cases, the 
models seem to underproduce the Si ll abundance with re- 
spect to Fe in comparison with observations, but since the 
yields used in the model are uncertain by a factor 2, at best, 
it is unlikely that the absolute value of the abundance ratio 
is known with an uncertainty better than 0.6 dex. 

In an alternative approach, Calura, Matteucci & Vladilo 

(2002) have produced models for a variety of galaxy types. In 
Figure 8, the light-coloured solid (dashed) curves are for spi- 
rals (irregulars). Indeed, in the chemical evolution models. 
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Figure 8. [Si/Fc] ratio as a function of mctallicity. [Si/Fe] is similar in both DLAs and sub-DLAs. The dark-coloured soUd (dashed) 
curves are contours for sub-DLAs (DLAs) column densities derived from models of disc galaxies (LIou, Boissier & Prantzos 2001). The 
light-coloured solid (dashed) curves are for spirals at 20 kpc (irregulars) representative of sub-DLAs (DLAs) in the models of (Calura, 
Matteucci & Vladilo 2003). 



irregulars are found to reproduce well the observed prop- 
erties of DLAs. In these models, there is no artificial bias 
applied to mimic the effect of dust on the observations, nor 
was it necessary to add an extra HI column density cut-off 
to match the quasar absorbers' definition. In the case of spi- 
ral galaxies, in particular, the models presented in Calura, 
Matteucci & Vladilo (2002) never reach the sub-DLA col- 
umn density since star formation is set to occur only above 
a density threshold, as suggested by Chiappini et al. (1997). 
In order to reproduce the properties of sub-DLAs, it has 
been necessary to eliminate this threshold and extend the 
computations towards larger radii (20 kpc). For [Si/Fe], the 
two types of models show little difference. In addition, the 
interpretation of this elemental ratio might be complicated 
by the effect of dust. 

[O/Fe]: O I is not often measured in DLAs, since the 
most predominant line, O I 1302, is often saturated. In con- 
trast, in sub-DLAs where the column densities are consid- 
erably smaller, this line is well suited for the abundance de- 
termination of oxygen. This element, which is not affected 
by dust depletion, is also a typical product of Type II su- 
pernovae. The prospect of measuring O I in large samples 
of sub-DLAs means that we will have a good indicator of a 
over Fe-peak element ratios. 

Figure 9 shows the evolution of the [O/Fe] ratio as a 
function of metallicity. All these ratios are recomputed as- 
suming [0/H]q — —3.26 (Holweger 2001). In most cases. 



this ratio is is sub-solar at high metallicities ([Fe/H]>-1.6). 
However, [O/Fe] is still greater than solar at low metallici- 
ties. Similarly, Galactic metal-poor stars show an enhance- 
ment between 0.35 and 1 (Goswami & Prantzos 2000 and 
references therein). However, in the case of high [O/Fe] ratio, 
one cannot decipher whether dust depletion or nucleosyn- 
thesis produces such high values. In the case of low [O /Fe] , 
however, we know that the intrinsic values cannot be much 
smaller than the values observed. 

Models from Hon, Boissier & Prantzos (2001) seem to 
underproduce the corresponding ratios with respect to the 
observations for both DLAs and sub-DLAs. In contrast, 
models from Calura, Matteucci & Vladilo (2002) seem able 
to reproduce the observed sub-solar values of [O /Fe] (see also 
Chiappini, Romano & Matteucci 2002). At present, however, 
the observed sample is still too small to allow firm conclu- 
sions to be drawn. 

[C/Fe]: For the same reasons as O I, the first compre- 
hensive set of measurements of C ll has been made in high 
HI absorbers. [C/Fe] is approximately solar in the majority 
of the Galactic stars. Previous studies of DLAs only allowed 
for lower limit determinations (see Figure 10), but measure- 
ments of C II in sub-DLAs provide a tracer of Fe-peak el- 
ements regardless of dust depletion. In this figure we use 
[C/H]q = —3.41 (Holweger 2001) as an estimate of the so- 
lar value. Some of the [C/Fe] ratios in sub-DLAs are clearly 
oversolar, at the level of the overabundance observed in the 
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Figure 9. [O/Fc] ratio as a function of metallicity. At low metallicity the ratio are greater than solar. Sub-DLAs generally have low 
[O/Fe], some of which are sub-solar. The dark-coloured solid (dashed) curves are contours for sub-DLAs (DLAs) column densities 
derived from models of disc galaxies (Hon, Boissier & Prantzos 2001). The light-coloured solid (dashed) curves are for spirals (irregulars) 
representative of sub-DLAs (DLAs) in the models of (Calura, Matteucci & Vladilo 2003). 



two tentative [C/Fe] measurements in the DLAs. Since C is 
a mildly refracted element in the interstellar medium, these 
measurements provide our best indicator of the dust con- 
tent of sub-DLAs. Indeed, although C and Fe do not share 
the same nucleosynthesis origin, in the Milky Way stars, C 
and Fe show identical behaviours, i.e. [C/Fe]=0 (Goswami 
& Prantzos 2000). Therefore the deviation from zero ob- 
served in sub-DLAs are suggestive of either the presence of 
dust in these systems or a different behaviour of C rela- 
tive to Fe relative to the Milky Way stars. Both types of 
chemical evolution models predict abundance ratios in the 
intermediate range between DLAs and sub-DLAs. Never- 
theless, the models of Hou, Boissier & Prantzos (2001) do 
not include yields from intermediate-mass stars. This will 
affect the [C/Fe] abundance ratio. The models of Calura, 
Matteucci & Vladilo (2002) predict very little difference be- 
tween DLAs and sub-DLAs, as is observed. But again, given 
the limited data set, these conclusions are tentative. 

[Al/Fe]: The metallicity evolution of the [Al/Fe] ratio 
is presented in Figure 11, using [A1/H]q = —5.51 (Grevesse 
& Sauval 1998). In DLAs, this ratio follows a trend simi- 
lar to the one observed in metal-poor stars (Prochaska & 
Wolfe 2002). In the sub-DLA case, this ratio is strongly af- 
fected by ionization corrections which increase the real abun- 
dance of Al with respect to that measured solely from Al ll. 
The [Al/Fe] ratio will consequently rise by a factor of 0.3 
dex in the worst case (see Paper I). Again sub-DLAs are 



characterised by a mean [Al/Fe] elemental ratio larger than 
DLAs. In this case, we note that models of Hou, Boissier 
& Prantzos (2001) do not match the values of high-redshift 
DLAs nor the observations in the Milky Way. It should be 
noted that the nucleosynthesis situation of this element is 
unclear (Goswami & Prantzos 2000) but it could also be 
that the intermediate-mass stars not included in the models 
contribute to the true [Al/Fe] abundance ratio. 



6 CONCLUSION 

Peroux et al. (2003) have extended the DLA definition to a 
new class of absorbers: the sub-DLAs with column density 
10^^ < N(III) < 2 xl0^° atoms cm"^ These systems are 
believed to contain a large fraction of the neutral gas mass 
in the Universe, especially a,t z > 3.5. Based on these consid- 
erations, we have constructed and fully analysed a sample of 
12 sub-DLAs (see Paper I). In the present paper, we analyse 
several of the properties of these absorbers in conjunction 
with DLAs from the literature. Our main findings can be 
summarised as follows: 

• Our sample of sub-DLAs can be used to observationally 
determine for the first time the shape of the column density 
distribution, f{N), down to N(HI) = 10^^ cm"^. The resuhs 
are in good agreement with the predictions from Peroux 
et al. (2003). An evolutionary study of f{N) is not possible 
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Figure 10. [C/Fe] ratio as a function of metallicity. Sub-DLAs are markedly more over-solar than DLAs. The dark-coloured solid 
(dashed) curves are contours for sub-DLAs (DLAs) column densities derived from models of disc galaxies (Hou, Boissier & Prantzos 
2001). The light-coloured solid (dashed) curves are for spirals (irregulars) representative of sub-DLAs (DLAs) in the models of (Calura, 
Matteucci & Vladilo 2003). 



with the current sample illustrating the need for a sample 
of sub-DLAs at high-redshifts. 

• An analysis of the clustering of sub-DLAs in our sample 
shows that these are more clustered than expected from a 
random distribution, although the number statistic is too 
small to draw firm conclusions. 

• We measure the kinematic properties of the sub-DLAs 
from our sample together with a statistically significant 
number of DLAs from the literature. We compare low- and 
high-ionization transition widths and find that the sub-DLA 
properties roughly span the parameter space of the DLAs, 
thus challenging multicomponent semi-analytical models, 
such as the one presented by Mailer et al. (2002). 

• The metallicity of absorbers as traced by [Fe/H] shows a 
slightly more pronounced slope for sub-DLAs {a — —0.40 ± 
0.22) than for DLAs {a = -0.18±0.12). In addition, the HL 
weighted mean metallicity is computed for various sub-sets 
of quasar absorbers. The evolution of [(Fc/Hdla)] might be 
stronger for sub-DLAs than for DLAs, and absorbers with 
N(HI) > lO^^" atoms cm ^ appear to be the less evolved, 
suggesting these objects in particular should be studied in 
detail. Observational evidence supports the hypothesis that 
this different behaviour is not due to the hidden effect of 
dust. We therefore propose that sub-DLAs might be associ- 
ated with a class of galaxies which better traces the overall 
chemical evolution of the Universe. 

• A study of the metallicity evolution with metal line 



profile ionization width might show hints of a correlation, 
whereby higher [Fe/H] ratios are associated with systems 
with larger widths. However the statistical significance of 
this result is low and a larger sample of observations is re- 
quired to unambiguously address this issue. 

• Abundance ratios for [Si/Fe], [O/Fe], [C/Fe] and [Al/Fe] 
were determined and compared with two different sets of 
models of the chemical evolution of galaxies. Overall, these 
appear to resemble abundance ratios observed in DLAs. The 
question thus remains open whether the two classes of ob- 
jects have similar chemical evolution histories or whether 
there are objects with different ages and star formation his- 
tories, but at a similar stage of chemical evolution. The first 
comprehensive sets of measurements of O I and C ll in high 
H I column density systems are given. Indeed, these are well- 
defined in sub-DLAs while they are almost always saturated 
in DLAs. These elements, unaffected by dust depletion, pro- 
vide direct indicators of the abundances in these systems. 



These various issues illustrate the importance of fur- 
ther studies of sub-DLAs to interpret the overall chemical 
evolution of neutral matter with redshift, as well as to ex- 
tend the analysis of quasar absorbers properties to a lower 
column density range. We are currently undertaking an ob- 
servational program aimed at studying sub-DLAs at high- 
redshifts. 
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Figure 11. [Al/Fc] ratio as a function of metallicity. Again, sub-DLAs are remarlcably more over-solar tlian DLAs. The solid (dashed) 
curves are contours for sub-DLAs (DLAs) column densities derived from models of disc galaxies (Hou, Boissier & Prantzos 2001). 
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